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ABSTRACT 



Context. Fossil galaxy groups are speculated to be old and highly evolved systems of galaxies that formed early in the universe and 
had enough time to deplete their L* galaxies through successive mergers of member galaxies, building up one massive central ellipti- 
cal, but retaining the group X-ray halo. 

Aims. Considering that fossils are the remnants of mergers in ordinary groups, the merger history of the progenitor group is expected 
to be imprinted in the fossil central galaxy (FCG). We present for the first time radial gradients of single-stellar population (SSP) ages 
and metallicites in a sample of FCGs to constrain their formation scenario. We also measure line-strength gradients for the strongest 
absorption features in these galaxies. 

Methods. We took deep spectra with the long-slit spectrograph ISIS at the William Herschel Telescope (WHT) for six FCGs. The ob- 
tained spectra are fit with Pegase HR SSP models within the full-spectrum fitting package ULySS yielding SSP ages and metallicities 
of the stellar populations. We measure radial gradients of SSP ages and metallicities along the major axes. Lick indices are measured 
for the strongest absorption features to determine line-strength gradients and compare with the full-spectrum fitting results. 
Results. Our sample comprises some of the most massive galaxies in the universe exhibiting an average central velocity dispersion of 
o"o = 271 ±28 km s~' . Metallicity gradients are throughout negative with comparatively flat slopes of V^e/H] = -0.19±0.08 while age 
gradients are found to be insignificant (V^ge = 0.00 ± 0.05). All FCGs lie on the fundamental plane, suggesting that they are virialised 
systems. We find that gradient strengths and central metallicities are similar to those found in cluster ellipticals of similar mass. 
Conclusions. The comparatively flat metallicity gradients with respect to those predicted by monolithic collapse (V z = -0.5) suggest 
that fossils are indeed the result of multiple major mergers. Hence we conclude that fossils are not 'failed groups' that formed with a 
top heavy luminosity function. The low scatter of gradient slopes suggests a similar merging history for all galaxies in our sample. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: stellar content - Galaxies: evolution - Techniques: spectroscopic 
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1. Introduction 

Fossil galaxy groups are observationally defined as extended X- 
ray sources with X-ray luminosities above Lx. boi ~ 

0.5 x 10 4i 

/i 5( , erg s _1 and a single bright elliptical dominating the opti- 
cal, all other galaxies within one-half virial radius bei ng at least 
2 magnitudes fainter in the R band (| Jones et al.|[2003h . Earlier 
numerical simulations of lBarnes1 (ll989) predicted such observa- 
tional characteristics for old and evolved groups of galaxies that 
had enough time to deplete their L* galaxies through successive 
mergers of member galaxies, building up one massive central el- 
liptical, but retaining the group-sized X-ray halo. Contrary to this 
hierarchical merging scenario it was also suggested that fossils 
could simply be 'failed groups' that for med with a top-heavy lu- 
mino sity function absent of L* galaxies dMulchaev & Zabludoffl 
[1991 . 

The first fossil group was identified by Ponma n et al. (1994). 
Since then, several other o bjects have been assigned to this 
class in the last few years dvikhlinin et al.1 Il999t iJones et alJ 
20031: lYoshioka et all 120041: 1 Sun et al.ll2004t IKhosroshahi et at] 
2004[ lUlmeretal l l2005USchkmer et all boiOt iPierini et al.l 
201 U iMiller et all 120121) . Dones et all (120031) estimated fossils 
to be quite frequent, constituting 8-20% of all systems with 
comparable X-ray luminosity (> 10 42 erg s _1 ). The detection 
of these aggregates is challenging however, mainly because of 
their unremarkable appearance in the optical and the lack of 
high S IN X-ray data. Recent studies have shown that fossils 



are much more massive than expected for remnants of group- 
sized haloes, comparable to po o r clus t ers rather than groups 



2006; 



(MendesdeOliveiraetal. 2006, 2009; Cvmiano et al. 
iProctor et al.ll2.01 it lEigenthaler & ZeiUngeHl2012l) . Because of 
these inconsistent findings, mainly due to the small number of 
fossils studied so far, various formation scenarios have been 
proposed for these system s. Based on scaling relat ions of or- 
dinary groups and clusters, IKhosro shahi et all (l2007|) concluded 
that fo ssils had an early formation epoch, while IProctor et all 
stated that fossils might simply be dark clusters, being 
comparable to ordinary clusters in mass and X-ray luminos- 
ity, but otherwise possessing the ric hness and optical lum inosi- 
ties of poor groups. Simulations by Dariush et ail ( 2007() using 
the M illennium and Millennium gas simulations fSpringel et al I 
2005) revealed that at any given redshift, fossils have assem- 
bled a higher fraction of their final halo mass compared to non- 
fossils with similar X-ray luminosity, suggesting that fossils in- 
deed formed early with most of their hal o mass in place at higher 
redshifts. A similar result was found by iD'Onghia etail (120051) 
using N-body/hydrodynamical simulations, concluding that fos- 
sils have already assembled half of their final dark matter mass 
at z > 1, leaving sufficient time for L* galaxies to merge, re- 
sulting in the exceptional magnitude gap at z — whilst in 
non-fossils the final dark matter halo mass is assembled quite 
recently. Although fossil systems are believed to have collapsed 
early and have assembled most of their virial masses at higher 
redshifts, the central galaxies in fossil groups may have merged 
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later t han non-fossil bright central galaxies (Dfa z-Gimenez et all 
l2008h . 

Assuming that fossils originate from ordinary galaxy groups 
that evolved to compact groups which subsequently coalesced 
into one massive elliptical, the merger history of the progeni- 
tor group is expected to be imprinted in the fossil central galaxy 
(FCG) and any differences in the evolution compared to field or 
cluster ellipticals should be reflected in the overall characteris- 
tics of these galaxies. Since the space density of fossils is con- 
sidered to be as large as that of poor and rich clusters combined, 
fossils could also provide suitable environments for the forma- 
tion of luminous cD gal axies observed in cluster s. Following 
this possible connection, iKhosroshahi et al.l d2006l) investigated 
the morphology of seven FCGs with the conclusion that the 
isophotal shapes of FCGs are always discy in contrast to bright- 
est gro up (BGG) and cluster (BCG) galaxies. lLa Barbera et al.l 
(2009) found no differences in the structural parameters between 
a sample of 25 FCGS, selected from the SDSS, and field el- 
lipticals. The examination of stellar populations revealed simi- 
lar ages, metallicities, and a-enhancement compared to bright 
field ellipticals. In addition, elliptical galaxies in compact groups 
were found to be older and more metal poor than field ellipticals 
and fos sil groups, which had already been e stablished in earlier 
works (iProctor et all 120041 IdelaRosa et al.ll2007l) . It was thus 
noted that FCGs cannot be formed by dry mergers of ellipticals 
in compact groups. However, a wet merger with a gas-rich disk 
galaxy could also explain th e observed relations. More recently, 
Mendez-Abre u et aTT(l2012l) have studied the near-infrared pho- 
tometric and structural properties of a sample of 20 FCGs. 
These fossils closely follo w the observed tilted fundamental 
plane of normal ellipticals (Bern ardi et al.l l2003bl) . suggesting 
that they are dynamically relaxed syste ms that suffered dissi - 
pational mergers during their formation (Hopkins et al . 2008). 
However, a change in the slope of scaling relations involving 
stellar mass was interpreted by subsequent major dissipationless 
mergers, increasing size and mass of the galaxi es but leaving 
velocity dispersion unaltered (Bernardi et al.ll20TTb . Because of 
low S/N data these authors couldn't determine whether fossils 
are systematically boxy or discy. Since most of these previous 
studies have focused on the photometric and structural properties 
of FCGs, detailed information on their spatially resolved stellar 
population parameters is still lacking. 

This paper investigates age, metallicity and line-strength gra- 
dients of a sample of six FCGs to shed light on the formation and 
evolutio n of fossils. The galaxies were s e lected from t h e sam - 
ples of iMendes de Oliveira et al.l d2006l) . ISantos et al] (120071) . 
and lEigenthaler & Zeilingen (120091) . according to their redshift 
and high apparent brightness. At redshifts between 0.023 < 
z < 0.081, the systems lie at luminosity distances between 
100 < Dl < 370 Mpc. These values correspond to an age of the 
universe of 13. 1 back to 12.4 Gyr at th e reds hift of the observed 
galaxies. According to iDariush et al.l (120071) . the Millennium 
simulation predicts that fossils at these redshifts should have as- 
sembled at least 95% of their total dark matter mass. 

The paper is organised as follows. In Sect. 2 observations 
and data analysis techniques are described whilst Sect. 3 presents 
all obtained results. Sect. 4 discusses and compares our findings 
with similar data from the literature. Conclusions are given in 
Sect. 5. Magnitudes presented in this work are SDSSIII DR9 
model magnitudes. Throughout the paper, the standard ACDM 
cosmology with Q.m = 0.3, Qa = 0.7, and a Hubble constant of 
Hq = 70 km s _1 is used. 



2. Observations and data reduction 

Long-slit spectroscopy of a sample of six FCGs was performed 
with the long-slit spectrograph ISIS at the William Herschel 
Telescope (WHT) in visitor mode. The observations were car- 
ried out in 2008 on December 20 and in 2009 on April 28. ISIS 
is suitable for long-slit spectroscopy up to 4 arcmin slit length. 
Through the use of dichroic filters, the instrument allows simul- 
taneous observations in both a blue and a red arm. Both arms 
provide an array of 4096x2048 pixels with 13.5^m and 15^m 
pixel size, respectively. The CCDs are arranged so that the long 
axis is placed along the dispersion direction maximising spectral 
coverage. The spatial scale is 0.2 arcsec pixel 1 in the blue and 
0.22 arcsec pixeL 1 in the red yielding a maximum unvignetted 
slit-length of 3.7 arcmin. The spectrograph is designed such that 
a slit- width matching the typical seeing of the site (~1 arcsec) 
corresponds to approximately 4 pixels FWHM on the CCD for 
all gratings. To cover the most prominent absorption features H/3, 
Mg 2 , Mgb, Fe5270, Fe5335, Fe5406, the 6100 dichroic in com- 
bination with the R300B grating (0.86 A pixeL 1 ) in the blue and 
the R1200R grating (0.26 A pixeL 1 ) in the red was requested. 
If no emission lines were expected around Ha from SDSS spec- 
tra, the R316R grating (0.93 A pixeL 1 ) was used instead. A slit 
width of 1 arcsec was considered for all galaxies. 

Table Q] lists all ellipticals observed at the WHT. The total 
integration time of the targets was split into 3 or 4 single ex- 
posures of ~ 45 minutes each to correct for cosmic rays. The 
spectrograph slit was aligned to the major axis of each galaxy 
except for RX J0752.7+4557 and RX J0844.9+4258 where the 
slit was aligned so that nearby bright galaxies were also covered. 
The slit position angles are given in Table Q] 

Basic data reduction and calibration was carried out in IRAF 
based on bias, flat, and arc frames taken in the afternoon and at 
dawn. Arcs were also taken during the night after each move- 
ment of the telescope to account for possible detector shifts. 
General CCD reduction steps such as bias subtraction, trim- 
ming, and flat-fielding were carried out within the IRAF ccdred 
package to remove the instrument signature. Bad pixels were re- 
moved by interpolating over nearby pixel regions. Subsequent 
spectral reduction was performed in the long si it package. 
Small geometric distortions introduced by the camera optics 
were removed by tracing the spectrum on the detector and trans- 
forming the tilted coordinates to the orthogonal image axes. 
Spectra were wavelength calibrated by identifying ISIS CuNe 
and CuAr lines in the arc lamp spectra. Line positions were fit by 
a 3 rd order cubic spline to set the dispersion solution. Individual 
frames were average combined to improve the S/N and filter 
out cosmic rays. Sky subtraction was carried out with the IRAF 
background task by specifying two or three background win- 
dows on either side of the combined galaxy spectra. A 6 th order 
Chebyshev polynomial was then fit to these background regions 
in all columns along the dispersion direction and subtracted from 
the data. The fit to the sky was checked for all galaxies in steps 
of 100 pixels along the dispersion direction. We note that be- 
cause of the small spatial extent of the investigated galaxies with 
respect to the CCD chip, the Chebyshev polynomial was always 
flat or showed only a very minor linear trend along the galaxy 
which was always much smaller than the noise of the sky back- 
ground. The background windows were chosen so that a large 
area of the sky was covered along the slit. The comparatively 
high degree of the Chebyshev polynomial was chosen so that 
a good estimate for the sky was achieved along the whole slit. 
The polynomial was almost always found to be flat. The sky 
subtracted spectra were then checked individually for any re- 
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Table 1. Fossil ellipticals observed with ISIS. 



galaxy 




c 

"2000 


7 


r' 
[mag] 


exposure 
[seconds] 


PA 

n 


spectral 
blue arm 


range 

red arm 


FWHM 

[arc sec] 


NGC 1132 


02 52 51.8 


-01 16 28.8 


0.023 


12.19 


4 x 2700 


140 


3800- 


-6250 A 


6600-7200 A 


1.3 


RX J0752.7+4557" 


07 52 44.2 


+45 56 57.4 


0.052 


14.43 


3 x 3000 


76 


3800- 


-6250 A 


6600-7200 A 


1.3 


RX J0844.9+4258" 


08 44 56.6 


+42 58 35.7 


0.054 


13.93 


3 x 3000 


36 


3800- 


-6250 A 


6600-7200 A 


1.3 


RX Jl 152.6+0328* 


11 52 37.6 


+03 28 21.8 


0.081 


14.37 


3 x 2700 


23 


3800- 


-6250 A 


6300-8500 A 


0.7 


RX J 1 520.9+4840* 


15 20 52.3 


+48 39 38.6 


0.074 


13.68 


3 x 2600 


94 


3800- 


-6250 A 


6300-8500 A 


0.7 


RX J1548. 9+0851" 


15 48 55.9 


+08 50 44.4 


0.072 


13.20 


3 x 2600 


33 


3800- 


-6250 A 


6300-8500 A 


0.7 



Notes: 

8 from 
* from 



Santos etal.ld2007h 



Eigenthaler & Zeilinger (2009) 
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Fig. 1. Observed spectra in the ISIS blue arm smoothed to the Lick/IDS resolution (~ 9 A). The most prominent spectral features 
and Lick/IDS index passbands (grey) are highlighted. The spectra corerspond to the central pixel of each galaxy. 
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maining night sky emission lines. For the fainter emission lines, 
the subtraction yielded an excellent result. Only the strongest 
lines, [OI] at /15577A and the sodium doublet at /15890A showed 
minor residual patterns. ING spectrophotometric standard stars 
matching the airmass of the scientific targets were selected from 
the ING landscape tooQ. Standard s tar spectra we re extracted by 
means of the algorithm described in lHornel(ll986l) to retrieve the 
maximum S/N without biasing the resulting fluxes. To minimise 
light loss due to differential atmospheric refraction, all stars were 
observed with the slit oriented along parallactic angle. Flux cali- 
bration was performed by comparing the observed spectrophoto- 
metric standard star count rates with the corresponding absolute 
flux tables in La Palma Technical Notes 65 and 100. To account 
for atmospheric extinction, the data from La Palma Technical 
Note 3 1 was applied. Due to the long exposure times, effective 
airmasses ha ve been calculat ed for each individual galaxy via the 
formula from Stetson (1989). The resulting sensitivity functions 
were finally applied to the galaxy spectra, removing both atmo- 
spheric extinction and the sensitivity characteristics of the CCD. 
To transform observed absorption-line features into the Lick IDS 
system, Lick standard stars were selected from the homepage 
of Dr. Guy WortheyQ. To match the stellar populations of the 
observed galaxies, G and K giants were considered for that pur- 
pose. FigureQ~]shows the final reduced spectra of the central pixel 
of each galaxy in the blue arm smoothed to the resolution of the 
Lick/IDS spectrograph. 

3. Data analysis and results 

A reliable determination of ages and metallicities in ellipticals is 
difficult, since both quantities alter spectra and colours of these 
galaxies in a similar way. An age change of +30% is cancelled 
by a simultaneous metallicity change of -20%. This is widely 
known as the age-metallicity degeneracy, or 3/2 degeneracy 
(IWorthevl fl999). Individual spectral features sensitive to either 
age or metallicity can be used to break this degeneracy, however 
dBurstein et al.ll 19841; IWorthev etal]|1994t IWorthev & Ottavianil 
1997). These well-calibrated features allow predictions for the 
absorption-line strengths of whole sing le stellar population s 
(SSPs) of a given age and metallicity (Worth ev et al.l [l994). 
When focusing only on specific features in the spectrum of an 
old stellar population, one obviously misses the wealth of in- 
formation embedded in all the remaining pixels along the dis- 
persion direction, however. To account for this additional infor- 
mation besides the strongest absorption features, full spectrum 
fitting techniques have been developed. Because of the very low 
S/N in absorption lines at larger radii of the observed ellipticals, 
ages and metallicities are solely based on full-spectrum fitting in 
this work. 

3.1. Full spectrum fitting with ULySS 

The open-source package ULyS SB (University of Lyon 
Spectroscopic Analysis Software; Kol eva et al.i r2009) was used 
to fit synthetic single stellar population (SSP) models directly 
to the observed galaxy spectra yielding SSP equivalent ages 
and metallicities of the stellar populations of the investigated 
galaxies. ULySS make s use of the PegaseHR models from 
iLe Borgne et al.l d2004|) based on the E lodie 3.1 library of stel- 
lar spectra dPrugniel & Soubiranll200Tl) . Since the synthesis of 

1 http : //catserver . ing . iac . es/landscape/ 

2 http : // astro . wsu . edu/wor they/html/ system . html 

3 |http : // ulyss . uni v- lyonl . fr/| 



a stellar population requires a stellar spectrum at any point in 
the parameter space (T e g, log g and [Fe/H]), the library in- 
volves an interpolator for that purpose. All Pegase HR models 
are computed assuming a Salpeter IMF and Padova 1994 evo- 
lutionary tracks providing synthetic SSPs with ages between 
1-20000 Myr and metallicities between -2.3 and 0.69 dex. A 
multiplicative polynomial is used to adjust the overall spectral 
shape to the SSP model. The order of the polynomial can be 
freely chosen but shouldn't be taken too low resulting mostly 
in a se vere mismatch of th e investigated spectrum and the SSP 
model. iKoleva et alj d2009h have shown that the model fits con- 
verge to stable values between 10 < n < 15 for the spectra of F, 
G, K and O stars and n > 35 for an M star and becoming unsta- 
ble again above n = 150. Hence, in order to achieve reasonable, 
stable values an order of 40 was chosen except where otherwise 
stated. 

Figure |2]presents the variation of the fit results with increas- 
ing order of n, specifically the deviations to values at n — 40. 
Shaded areas correspond to deviations of 0.05 dex in [Fe/H] 
and 0.1 dex in age. Except for RX J1548. 9+0851, values stay 
within 0.05 dex in [Fe/H] from 40 < n < 80, matching the av- 
erage error of the measured values. For orders above n = 80, 
values start to diverge again since the polynomial is trying to 
over-fit the spectra. Values of age show a larger spread com- 
pared to metallicity. The variations of age and metallicity with 
increasing n are likely due to the limited S/N of the observed 
spectra. The ULySS parameter \CLEAN was considered for the 
fitting to exclude possible outliers in the galaxy spectrum, re- 
sulting from any remaining minor night sky emission. The clip- 
ping algorithm reduces the sensitivity of a fit to such undesir- 
able features by reiterating the fit after outliers in the residu- 
als have been identified and maskecQ. The reliability of the per- 
formed SSP fits was tested through the use of x 2 and conver- 
gence maps in the age-metallicity parameter space. To deter- 
mine integrated SSP equivalent ages and metallicities of each 
galaxy, one-dimensional spectra were created by summing up 
all the signal along the spectrograph slit. 200 Monte-Carlo sim- 
ulations were computed for every galaxy, repeating a fit suc- 
cessively with random Gaussian nois^l applied to the spectrum 
in each step. The resulting point distributions were then used 
to calculate average SSP ages and metallicities and the associ- 
ated standard deviations. Outliers clearly detached from the main 
point distributions were excluded for these measurements. Only 
blue arm spectra were considered because of the lack of promi- 
nent absorption features in the red channel. Convergence maps 
confirmed the independence of the obtained results from varying 
initial guesses for age and metallicity. Figure [5] shows the final 
point distributions of age and metallicity resulting from the ap- 
plied Monte-Carlo simulations of the one-dimensional SSP fits. 
While the galaxies spread from 8 to 17 Gyr in SSP model age, 
SSP metallicities show consistent subsolar values around -0.15 
dex. Almost all Monte-Carlo simulations follow a similar trend, 
revealing a much higher uncertainty in age than in metallicity. 

3.2. Age- and metallicity gradients 

To investigate gradients of age and metallicity along the galaxy 
major axis, the observed long-slit spectra were chopped into bins 



4 see |http: //ulyss .univ- lyonl .fr/uly. html#ULY_FIT for a 
detailed decription. 

5 The dimension of the added noise was based on a user-defined 
signal-to-noise ratio provided for each galaxy. The S/N ratios were mea- 
sured with IRAF splot at a rest-frame wavelength of ~ 5 140 A. 
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Fig. 2. Variation of SSP age and metallicity with increasing order n of the multiplicative polynomial for all observed ellipticals. The 
deviations to values at n — 40 are displayed. Deviations corresponding to 0.05 dex in [Fe/H] and 0. 1 dex in age are shaded. 



along the spatial axis. The individual spectra were then fit sepa- 
rately with ULySS. Each spectrum was binned along the spatial 
direction to achieve a higher S/N. The seeing was estimated to 
1.3 and 0.7 arcsec for the two nights corresponding to about 7 
and 4 pixels, respectively. Central galaxy spectra were binned 
only by this small number of pixels, matching the present seeing 
while 30 - 40 pixel rows were binned in the outermost spec- 
tra to achieve a reasonable S/N along the whole slit. Again, the 
S/N was estimated for each spectrum with splot at a rest-frame 
wavelength of ~ 5 140 A. The S/N varies from ~ 60 - 25 for 
the central bins to ~ 19 - 10 for the outermost bins. Figures [3] 
and |4] show spectra of the central and outermost bins together 
with the corresponding SSP model fits and residuals. The esti- 
mated noise strongly affects the resulting Monte-Carlo simula- 
tions. While the central bins showed well-defined point distribu- 
tions, the outer ones yielded much more diffuse results. In addi- 
tion, there is a clear tendency that the fits in the outer bins tend 
to converge to the model limit at 19.99 Gyr. These data points 
were all neglected for the age and metallicity measurements. In 
each bin, the resulting ages and metallicities were then computed 
as the mean of the 200 Monte Carlo simulations. At a given 
radius, ages and metallicites were averaged from both sides of 
the galaxy centre. The corresponding error bars were estimated 
from the associated standard deviations. Assuming power-law 
gradients in age and metallicity, linear relations between log age, 
[Fe/H], and log r/r e ff are expected. Hence, linear relations of the 
form 



log (age[r/r eff ]) = log (age eff ) + V age log (r/r eff ) (1) 
[Fe/H] [r/reff] = [Fe/H] eff + V [Fe/H] log (r/r eff ) (2) 

have been fit to the SSP ages and metallicities obtained with 
ULySS, where V age and Vrp e /H] are the corresponding gradients 
in age and metallicity. The fits are presented in Fig. [6] while 
the corresponding numbers and input parameters for the fits are 
given in Table [2] 

All metallicity gradients exhibit negative slopes between 
-0.11 < V|Fe/H] ^ -0.31 while gradients in age show both 
positive and negative slopes between -0.05 < V[p e /H] ^ 0.06. 
Gradient errors were determined directly from the linear least 
squares fits. To quantify the goodness of the fits, the coefficient 
of determination R 2 is listed in Table [2] for all fits. A value of 
R 2 — 1 corresponds to a perfect linear relation while R 2 — in- 
dicates no linear relation at all. The observed metallicity profiles 
are much better fit by the utilized power-laws than the profiles in 
age. 
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Fig. 5. ULySS fit results for the integrated stellar populations 
of all observed FCGs. The point distributions result from 200 
Monte-Carlo simulations. Adopted final values and errors for 
age and metallicity are the averages and standard deviations 
of the shown point distributions indicated by the bold crosses. 
Obvious outliers have been excluded for the measurements. 



3.3. Line-strength gradients 

Besides SSP ages and metallicities, the strongest absorption 
lines have been measured to investigate line-stre ngth gradients. 
For th at purpo se, Lick/IDS indice s as de fined in IWorthev et al.l 
(1 19941) and IWorthev & Ottavianil (1 19971) have been analysed. 
Due to the limited S/N at larger galactic radii, only the most 
prominent indices (HJ3, Mg 2 , Mgb, Fe5270, Fe5335, Fe5406) 
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Fig. 3. ULySS SSP fits to the central bin of the observed FCGs. The galaxy spectra are shown in black while the best-fit SSP models 
are shown in blue. Fluxes are given in 10~ 17 erg s _I cirT 2 A -1 . S/N ratios of the investigated spectra are also presented. 



have been considered. Figure Q] highlights the redshifted index be degraded in resolution first. For that purpose, the FWHM 
bands in all six ellipticals. The presented spectra have been of emission lines was measured in the arc comparison spec- 
smoothed to the Lick/IDS resolution. Other prominent lines tra yielding an estimate of the actual spectral resolution. The 
are also marked. Since the spectral resolution of the Lick/IDS spectra were then convolved with a Gaussian kernel of disper- 
system amounts to only ~ 9 A, the observed spectra had to sion cr smooth = [(FWHM 2 ick - FWHM 2 bs )/81n2] I/2 to match 
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the Lick/IDS system. The broadened spectra we re subsequentl y timated at rest-frame wavelength of ~ 5140 A and 200 Monte 
analysed with the open-source program index j^ dCardiell2 007). Carlo simulations were considered. Final indices and errors were 
To estimate errors of the measured indices, again a fixed S/N es- then computed as the mean of the 200 Monte Carlo simula- 
tions. To transform the observed indices to the Lick/IDS system, 

— ; Lick standards have been observed. Figure |7]shows the compar- 

. version 3.0 - 24.08.2005 i son between the measured indices and the original Lick data 

http : //www . ucm . es/inf o/Astrof / software/ indexf/ indexf . html 
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Table 2. SSP equivalent central ages and metallicites, central velocity dispersions and gradients in age and metallicity obtained 
with ULySS. The corresponding fit parameters are also given. 



galaxy 


MD" 


LMIN 


-LMAX 6 


age [Gyr] 


[Fe/H] 


cr [km s ']'' 


V [Fe/H] d 




Vagi 


d 


R le 


NGC 1132 


50 


3800- 


6100 A 


15.1 + 1.0 


-0.11+0.03 


227+ 3 


-0.14 


± 0.02 


0.68 


-0.03 ± 


0.02 


0.13 


RX J0752.7+4557 


40 


3800- 


•5800 A 


9.7+1.8 


-0.01+0.06 


296+ 18 


-0.28 


± 0.17 


0.57 


0.00 ± 


0.08 


0.00 


RX J0844.9+4258 


60 


3800- 


•5800 A 


17.3 + 1.3 


-0.05 + 0.03 


284+ 4 


-0.12 


± 0.02 


0.87 


-0.04 ± 


0.02 


0.48 


RX Jl 152.6+0328 


40 


3800- 


•5800 A 


12.4+1.3 


-0.07 ± 0.04 


246+ 5 


-0.31 


± 0.14 


0.62 


0.06 ± 


0.03 


0.59 


RX J 1520.9+4840 


40 


3800- 


•5800 A 


11.1+2.7 


0.02 ± 0.07 


279+ 9 


-0.20 


± 0.06 


0.71 


-0.05 ± 


0.08 


0.10 


RX J1548. 9+0851 


40 


3800- 


-5800 A 


11.5 + 2.1 


-0.14 + 0.06 


292 ± 10 


-0.11 


± 0.09 


0.30 


0.04 ± 


0.02 


0.55 



Notes: 

" Order of polynomial degree. 

Wavelength range used for the SSP fit. 
c Aperture corrected central velocity dispersion. 
d Age- and metallicity gradients as defined in equations[T]and|2] 

e Coefficient of determination: R 2 = 1 - perfect linear relation; R 2 = - no linear relation. 



1.2 

1.0 
1.2 

1.0 

0.8 
1.4 



1.2 

O 

S 5 1.2 



1.0 
1.2 
1.0 



i i 



NGC1132 



_i I i i i i L_ 



-p— i— i— i— i— p 

RXJ0752.7+4557 



1— r 



i - 



~ I 1 1 T" 



1 ■ 1 l_ 



-i — I — i — r— 



-l — i — i — i - 



■ ■ ■ i_ 



~ I — i — i — i - 



1 — i — i — i — i — | — i — r- 

RXJ0844.9+4258' 



-I 1 1 1 - 



1 ■ 1 I— 



-i — i — i — r- 



■ ■ ■ i_ 



- i — i — i — i — | — i — r— 

RXJ1 152.6+0328 



_i i i i_ 



- 1 — i — i — i 



RX J1 520.9+4840 - 



1.2 - 

1.0 - 1 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

RXJ1 548.9+0851 - 



i J f 



f i 



J I I I I I I I I I I I I I I I I I I I l_ 



-2.0 



-1.5 



-1.0 
log (r/r ) 



-0.5 



0.0 



00 ^ L— * 



-0.5 
0.0 

-0.5 r> 



0.0 
-0.5 
0.0 
-0.5 

0.0 
-0.5 
0.0 

-0.5 - 



i i i i i i i i i i i i 



i 

NGC1132 



A -*^-*-t-i.«iL 



_! I I I I I I I I I I 1_ 



- 1 1 1 1 1 1 1 1 1 1 - 



r-i— i— i— i— j— t— i— i— i— r— i— r- 
RXJ0752.7+4557 



>~~-L 



—I 1 1 1 1 1 1 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l~~ 

RXJ0844.9+4258 



I L 



-1 1 1 1 T" 







-i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

RXJ1 152.6+0328 J 



J I 1 L 



1 I 1 



■ ■ I 



J I I I I I I 1_ 



_l I I I I I l_ 



_l_ 



8 ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

RX J1 520.9+4840 



—i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i - 

RXJ1 548.9+0851 



j i i I i i i i I i i i i I i i i i I i_ 



-2.0 



-1.5 



-1.0 
log (r/r) 



-0.5 



0.0 



Fig. 6. Radial gradients of age and metallicity for the observed FCGs. Left panels refer to age, right ones present metallicity. Dashed 
lines show the linear relations fit to the observed data. The ordinates are the decimal logarithm of age and the metallicity in dex. 
Each datapoint is the result of 200 Monte Carlo simulations and was averaged from both sides of the galaxy centre. Error bars are 
the corresponding standard deviations. Metallicity gradients are throughout negative (-0.31 < V[p e /H] ^ -0.1 1) while age gradients 
are flat scattering around zero (-0.05 < V age < 0.06). Vertical lines indicate the seeing-dominated central parts. 



from IWorthev & Ottavianil (Q997) for the observed standards. 
Indices show only small differences compared to the original 
Lick data, proving the reliability of our measurements. To quan- 
tify the deviations to the Lick system, linear relations have been 
fit to the data. To account for Doppler broadening, indices were 
measured consecutively in the Lick standard star spectra broad- 
ened to 50-400 km s _1 in steps of 50 km s . Assuming that 
(index !Jcr ) and (index ,jo) are the absorption line-strengths of 
index i and standard star j measured at velocity dispersion a 



and 0, then % ; )(T = n 1 Yj"=i [index, j i0 - /index, jo] is the average 
correction factor estimated from all observed standard stars for 
index i. Thus, an index in a galaxy with velocity dispersion <x 
was corrected by means of index, o = index ( - ■ 9?7o- Not all in- 
dices show the same dependency with <x. While iron lines show 
a rather steep decline in index strengths, H/3 and Mg2 show com- 
paratively small dependence. Correction functions 9?, jCr were es- 
timated by averaging individual correction functions 7? y jCr (ex- 
cluding obvious outliers) and fitting third order polynomials to 
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12 3 4 

observed 



1 2 3 

observed 



Fig. 7. Comparison bet ween the measured lick ind ices and the 
original lick data from lWorthev & Ottavianil (1 1997b for the ob- 
served lick standards. Horizontal axes present the measured in- 
dices from this work while vertical ones show the corresponding 
values from the literature. Dashed lines indicate unity while solid 
ones present linear least squares fits to the data. 



the averaged data. Ultimately, all measured indices were cor- 
rected for velocity dispersion first and subsequently transformed 
to the published lick data. Gradients V/ of individual indices / 
were measured by linear least squares fits of the form 

/(r/ieff) = 7(r eff ) + V7 log [r/r eff ] (3) 

so that VI = A//Alog (r/r e ff). As for ages and metallicites, 
indices at a certain radius were averaged from both sides of the 
galaxy centre. Obvious outliers at larger radii deviant from the 
overall gradient trend were excluded for the fit. Figure |8]shows 
the resulting gradient slopes V7 with respect to the central veloc- 
ity dispersions of each galaxy. Dashed lines show the avereage 
line-strength gradient for all galaxies while short-dashed lines 
highlight flat gradients for comparison. In the case of Mgb, the 
very steep gradient for RX J1548. 9+0851 has been excluded for 
the estimation of the average gr adient slope. We comp a re our 
ind ex gradients with th ose from jKobavashi & Arimoto (1999) 
and iRawle et alj (|2010) and find similar gradient slopes. Only 
for Mgb and Fe5335 we find gradients that are much steeper 
than in our comparison sample. Table [3] lists the line-strength 
gradients for each index and also gives the coefficient of deter- 
mination for each fit. All metal sensitive lines (Fe and Mg) show 
negative gradients confirming the negative metallicity gradients 
derived from the SSP fits. Only Fe5406 in RX1548.9+0851 ex- 
hibits a flat gradient. The age sensitive H/3 index exhibits both 
negative and slightly positive gradients. The steep positive gra- 
dient for H/3 in RX J0752.7+4557 results from the fine being 
observed in emission. Hence, this slope is not shown in Fig. 
[8] The observed line-strength gradients follow the trend derived 
from the SSP fits, clearly indicating negative metallicity gradi- 
ents. To compare our index gradients with the metallicity gradi- 
ents from ULyS S, we consider index-metallicity relations from 
IWorthevi (1 19941) SSP mod els. In particular, we fit linear rela- 
tions to the lWorthevl(ll994 indices Mg 2 , Mgb, Fe5270, Fe5335, 
Fe5406 for metallicities -0.5 < [Fe/H] < 0.5 and galaxies with 
12Gyr. These relations are of the form [Fe/H] = p ■ index - d so 
that an index gradient VI corresponds to a metallicity gradient 



V[Fe/H] = p • VI. We derive values of p 



Ms!. 



5.48, p 



Mob 



0.44, 



PFe5270 = 0.63, /? Fe 5335 = 0.52, and p Fe 5406 = 0.79. Applying 
these conversion factors to the average line-strength gradient 



slopes illustrated in Fig. [8] we find average metallicity gradi- 
ents Of Vrpe/H] = -0.28 ± 0.05 (Mg 2 ), Vrpe/H] = -0.44 ± 0.22 
(Mgb), V [Fe/H] = -0.29 ± 0.1 1 (Fe5270), V [Fe / H ] = -0.29 + 0.24 
(Fe5335), V [Fe/H] = -0.26 + 0.18 (Fe5406) for the observed 
fossils. Except for the systematically lower value derived from 
Mgb, all gradient slopes show very similar values resulting in an 
average metallicity gradient of V[p e /H] = -0.3 1 +0.07. This result 
is consistent with the average gradient slope V[p e /H] = -0.19 + 
0.08 derived from full-spectrum fitting within the error bars. 
Individual gradient slopes are also consistent with the metallic- 
ity gradients derived from full-spectrum fitting within the corre- 
sponding error bars except for NGC1 132, showing a systemati- 
cally flatter gradient. We derive V[p e /H],NGCii32 = -0.26 + 0.05, 

V[Fe/H], RXJ0752 = -0.39 ± 0.15, V[ Fe /H], RXJ0844 = -0.22 ± 0.12, 
V[Fe/H],RXJ1152 = -0.28 ± 0.16, V [Fe /H], RXJ1520 = -0.41 + 0.16 

and V[Fe/H],RXJi548 = -0.18 + 0.12 from the averages of all line- 
strength gradients. However, we focus our results on the full- 
spectrum fitting technique since every pixel is considered for the 
determination of metallicity instead of only a few individual in- 
dices. 



3.4. Velocity dispersions 

To estimate the dynamical mass of the observed ellipticals, line 
of sight velocity dispersions (LOSVDs) have been derived with 
ULySS. By fitting SSP models to the observed spectra, the pack- 
age also measures the broadening of each SSP model to match 
the observed spectrum, allowing for a proper determination of 
LOSVDs. ULySS determines the broadening relative to the SSP 
model dispersion cr m(M ] e i, amounting to 13 km s _! for the Pegase 
HR models. In addition, the instrumental dispersion cr instl of the 
spectrograph was taken into account by measuring FWHMs of 
emission lines in the arc lamp comparison spectra taken during 
the night in the same setup as the galaxies. Based on the av- 
erage of several lines a final value of <Xj nstr = 95 km s _1 was 
considered^. Final LOSVDs were ultimately computed via the 
relation: a 2 , — a 2 . + a 2 , , - a 2 , . The observed central 

phys ulyss model mstr 

velocity dispersions have been measured in the inner seeing - 
limited bin of each galaxy. This central aperture corresponds 
to different physical radii based on the distance of each galaxy. 
To account for this inconsistency, we applied an aperture cor- 
rection, transforming the observed dispersions to a system that 
is independent of distance. For that purpose, the prescription 
C"0,corr — 0"obs • [8(r apelture /r e )] 04 has been appli ed, were 

^"aperture 

is the central aperture and r e , the effective radius (J0rg ensen et al.l 
119951; IWegner et ail 119991) . Table |2] shows the resulting, aper- 
ture corrected central velocity dispersions cr con . The integrated 
velocity dispersions cr were used to correct the measured ab- 
sorption indices for Doppler broadening. The studied FCGs ex- 
hibit an average central velocity dispersion of 27 1 ± 28 km s _1 . 
Com paring this value with the LOSVDs of ~ 9000 ellipticals 
from lBernardi et aLl(l2003cb shows that the investigated galaxies 
are among the most-massive ellipticals in the universe. 



3.5. The fundamental plane 

Elliptical galaxies show a tight relation between their effec- 
tive radius r e , the mean effective surface brightness {p) e - 
-2.5 log (I) e , and the central vel ocity dispersion crp known 
as the fundamenta l plane (FP) dDiorgovski & Davisl 11987b 
iDressler et al.ll987l) : r e oc cr$ (I e ) B , with varying values of A and 
B by different authors (see lMagoulas et al1l2012t iBernardi et al.l 



FWHM = 2 V21n2 • a 
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Table 3. Line-strength gradients as defined in equation[3] 



galaxy 


VH/3 


R 2 " 


V Fe5270 


R 2a 


V Fe5335 


R 2 " 


VFe5406 




VMg 2 




VMgb 


R 1 " 


NGC 1132 


0.0 = 


= 0.1 


0.02 


-0.6 = 


= 0.2 


0.56 


-0.4 = 


= 0.2 


0.32 


-0.3 = 


= 0.1 


0.56 


-0.05 = 


tO.Ol 


0.78 


-0.6 = 


= 0.1 


0.61 


RX J0752.7+4557 


1.8 = 


= 0.1 


0.99 


-0.4 = 


= 0.8 


0.09 


-1.0 = 


= 0.3 


0.81 


-0.7 = 


= 0.2 


0.81 


-0.04 = 


= 0.02 


0.73 


-1.0 = 


= 0.2 


0.91 


RX J0844.9+4258 


-0.2 = 


fc0.2 


0.24 


-0.2 = 


= 0.1 


0.43 


-0.1 = 


= 0.4 


0.02 


-0.3 = 


= 0.1 


0.69 


-0.06 = 


tO.Ol 


0.91 


-0.7 = 


= 0.1 


0.95 


RX Jl 152.6+0328 


-0.4 = 


= 0.5 


0.19 


-0.7 = 


= 0.5 


0.42 


-0.1 = 


= 0.6 


0.01 


-0.2 = 


= 0.3 


0.16 


-0.06 = 


= 0.03 


0.72 


-0.9 = 


= 0.4 


0.72 


RX J1520.9+4840 


-0.3 = 


= 0.3 


0.24 


-0.6 = 


= 0.4 


0.30 


-1.2d 


= 0.4 


0.74 


-0.5 = 


= 0.3 


0.41 


-0.05 = 


= 0.02 


0.65 


-1.8 = 


= 1.1 


0.41 


RX J1548. 9+0851 


0.0 = 


= 0.2 


0.01 


-0.4 = 


= 0.2 


0.49 


-0.5 = 


= 0.3 


0.30 


0.0 = 


= 0.2 


0.00 


-0.04 = 


= 0.02 


0.64 


-4.9 = 


= 2.1 


0.57 



Notes: 

" Coefficient of determination: R 2 = 1 - perfect linear relation; R 2 = - no linear relation. 
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Fig. 8. Line-strength gradients versus c entral velocity dispersion. Bla ck tri angles show the o bserved FCGs while grey squares and 
triangles are index measurements from iKobavashi & Arimotol (1 1999I) and iRawle et alJ d2010l) . Dashed lines represent the average 
gradient for all observed fossils while short-dashed lines indicate zero gradients. The HB gradient is not shown for RX J0752. 7+4557 
since the line is seen in emission. The very steep Mgb gradient of RX J1548. 9+0851 has been excluded for the estimation of the 
average Mgb gradient. 



2003b). Since this fundamental plane spans about three orders of 
magnitude in luminosity, exhibiting very low residual scatter and 
little variation with environment, it suggests a uniform formation 
process of ellipticals. The FP can be reconstructed by means of 
the virial theorem which implies values of A = 2 and B — -1 
and a constant mass-to-light ratio along the plane. However, ob- 
served values of A — 1.49 and B - -0.75 from a sample of ~ 
9000 SDSS early-type galaxies bv IBernardi et alJ (l2003bl) indi- 
cate a tilted plane that deviates from the virial prediction. 

To study the location of the observed FCGs on the FP 
we modeled the SDSS III DR9 surface brightn ess profiles with 
GALFIT0 by fitting a Ide Vaucouleursl (1 19481) r 1/4 law to the 
galaxy surface brightness distribution. Nearby stars were used 
to estimate the PSF and correct for seeing. By applying a 
Ide Vaucouleursl (1 19481) r 1 ^ 4 law, the mean effective surface 
brightness (n) e can be computed as (fj) e — fi e - 1.393 (see 
Graham & Driver 2005). The measured surface brightness was 
corrected for (1 + zf cosmological dimming, extinction and 
.^-correction, taken from the SDSS. Effective radii were cor- 
rected for ellipticity so that r e = r e jQ b s yfbja, where b and a 
are the semi-minor and semi-major axes, respectively. As men- 
tioned in Sect. 13. 41 we applied an aperture correction to our cen- 
tral velocity dispersions. The same procedure was applied in 



IBernardi et al.l(l2003al) allowing us to directly compare our mea- 
surements with their work. 

Figure|9]sho ws the obser ved FCGs on the fundamental plane 
from Bern ardi et al.l d2003bl) . obtained from ~ 9000 early-type 
galaxies in the SDSS in the redshift range 0.01 < z < 0.3 
(solid line). The dotted lines present 3 cr deviations from the 
fit. All FCGs are located within these thresholds, following the 
SDSS fundamental plane. This result supports the idea that these 
objects are relaxed systems that had enough time to virialise. 
However, all galaxies are found below the FP, indicating that 
they are either more massive or more compact compared to the 
majority of ellipticals. Given that these galaxies are expected to 
have accumulated enough mass from L* galaxies of the progen- 
itor group, it is plausible that the observed fossils show system- 
atically higher masses compared to average ellipticals. 



4. Discussion 

Radial gradients of age, metallicity, and individual ab- 
sorption line-strengths have been exte nsively studied in 



the past for seve r al elliptical galaxies (|Davies et aL 1993; 



Fisher et al. 1995; Sanchez-Blazquez et al. 2006; Redaetal. 
[20071; lAnnibali et alJl2007t iBaes et alJl2007t IRawle et al.l l2010; 

lLoub ser & Sanchez-Blazquezl 120121 and references therein). 

http://users.obs.carnegiescience.edU/peng/work/galfit/gha:SagrhdMrits are in general well fit by power-laws and their 
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2.0 2.5 3.0 

log g + 0.20 «u> e - 20.09) 

Fig. 9. The observed FCGs on the fundamental plane. The solid 
line shows the best-fit FP from lBernardi et aD (|2003b) obtained 
from ~ 9000 early-type galaxies in the SDSS in the redshift 
range 0.01 < z < 0.3 while the dotted lines present the 3 <x 
deviations from the fit. Error bars for effective radii and surface 
brightness were determined with GALFIT while velocity disper- 
sion errors were estimated with ULySS. 



slopes provide useful constraints on the formation scenario of 
ellipticals. Following this approach, we have observed six fossil 
central ellipticals (FCGs) spectroscopically to study their stellar 
populations and measure age, metallicity, and line-strength gra- 
dients of the strongest absorption features to shed light on the 
formation and evolution of fossil systems. For that purpose we 
obtained medium-resolution long-slit spectra at the WHT and 
made use of the full-spectrum fitting package ULySS to fit SSP 
models to our data. 

Figure [10] shows central velocity dispersions, ages, and 
metallicities of our FCGs (red triangles) in com parison with 
the sa mple of early-type galaxies presented in iKoleva et al.1 
d201 lb . Black squares are E and SO galaxies in the Virgo and 
Fornax Clusters (see also ISpolaor et al.l l2009t iBedregal et al.l 
l2006UVazdekis etal. 112001 1) while g rey symbols sh ow early-type 
dwarfs from lDe Riicke et al1(l2005l) . The sample of Kole va et al.l 
(1201 ll) provides an ideal comparison for our measurements since 
it covers a broad range in galaxy mass and stellar populations 
are analysed as highlighted in this work, i.e. all galaxies are fit 
with Pegase HR models within ULySS. Exhibiting an average 
central velocity dispersion of <x = 271 + 28 km s _1 , our sam- 
ple comprises some of the most-massive ellipticals in the uni- 
verse, populating the high end mass-regime of early-type galax- 
ies. This is not surprising given the fact that all observed ellip- 
ticals are expected to be the remnants of the merger of multiple 
L* galaxies. The mass-metallicity relation is evident in Fig. \W[ 
more massive galaxies retaining their metals much easier than 
less-massive ones. We find similar central metallicities around 
[Fe/H]o = -0.06 + 0.06 for our fossils compared to the cluster 
ellipticals in Virgo and Fornax of similar mass. Central ages are 
old and exhibit an average value of 12.8 + 0.8 Gyr. Except for 
one galaxy, all observed ellipticals exhibit SSP ages older than 
10 Gyr indicating formation epochs at z > 2. RX J0752. 7+4557 
shows an intermediate age of 8.0 + 1.0 Gyr suggesting a more 
extended star formation history to z — 1 . This galaxy is also the 
only one in our sample with prominent nebular emission lines 
see Fig. [0 and i t was classified as AGN by multiple authors 
Hess et alJl20!2t IVeron-Cettv & Veronll2010l) . Metallicity gra- 



dients in our sample are found to be throughout negative with an 
average slope of V[p e /H] = -0.19 + 0.08 while age gradients are 
fiat, scattering around V age = 0.00±0.05. FiguresfTOlandfTTIshow 
that the observed gradients are similar to those of our control 
sample, suggesting that the physical processes involved in the 
formation of fossils do not differ systematically from those form- 
ing cluster ellipticals. The full-spectrum fitting technique ap- 
plied here is a powerful tool in breaking the age-metallicity de- 
generacy of old stellar populations. However, to confirm our SSP 
fitting results, we also measured radial gradients for individual 
line-strengths. Metal-sensitive aborption features such as Mg2, 
Mgb, Fe5270, Fe5335 all show negative radial gradients while 
the age-sensitive Hfi index shows both positive and negative gra- 
dient slopes, scattering around a flat gradient. We compare our 
ind ex gradients with t hose from [Kobayashi & Arimoto (1999) 
and lRawle et alJ (l2010h and find similar gradient slopes. Only for 
Mgb and Fe5335 we find gradients that are much steeper than in 
our comparison sample. From our index measurements we de- 
rive an average metallicity gradient of V[ Fe /H] = -0.31 + 0.07. 
This result is lower but consistent with the average gradient 
slope V[Fe/H] = -0.19 + 0.08 derived from full-spectrum fitting 
within the error bars. Individual gradient slopes are also con- 
sistent with the metallicity gradients derived from full-spectrum 
fitting within the corresponding error bars except for NGC1 1 32, 
showing a systematically flatter gradient. 

Two competing formation scenarios have been proposed 
for the formation of elliptical galaxies. Either ellipticals are 
forme d monolithical l y by dissipation al collapse of a pristine gas 
cloud dLarsonl["l974t ICarlbergill984l) or hierarchically through 
successive mergers of disk galaxies or many dwarf galaxies 
dKauffmann & Wh"itelll993l) . In the dissipative collapse model, 
gas gets continuosly enriched in the course of star-formation as 
it flows towards the galaxy centre, resulting in steep metallicity 
gradients. In this scenario gradient slopes have be en predicted t o 
be as steep as V? = A log Z/ A log r - -0.35 (Larson 1974), 
V z = -1.0 (iLarsonl 11975k . and V z = -0.5 dCarlberg| [i984). 
Indeed, negative metallicity gradients, observed as radial gra- 
dients of absorption line-strengths, are a common feature in el- 
lipticals, however exhibiting typically shallower slopes around 
Vz = -0.3. Monolithic collapse also predicts a strong correla- 
tion between galaxy mass and metallicity gradients since deeper 
potential wells retain more metals. 

In the framework of hierarchi cal ma s s asse mbly, ellipticals 
are the results of galaxy mergers. White (1980) has shown that 
mergers flatten and erase the metallicity gra dients of the progen- 
itor g alaxies due to orbital mixing of stars. Mih os & Hernquisn 
(1 1994 modelled the merger of disk galaxies and concluded that 
pure stellar disks do not reproduce the power-law gradients ob- 
served in ellip ticals unless a metal-rich starburst is involved. 
KobavashJ d2004l) studied the evolution of metallicity gradients 
in the framework of ACDM cosmology, in which the formation 
and evolution of ellipticals is expected to be driven by merg- 
ers. Allowing for various merging histories for ellipticals, rang- 
ing from quasi-monol ithic al collapse to multiple major merg- 
ers, Kobayashl d2004l) found no correlation between gradient 
slopes and galaxy mass and stated that the variety of gradients 
originates in the different merging histories of ellipticals in the 
sense that galaxies that form monolithically have steeper gradi- 
ents than galaxies resulting from major mergers. The change in 
gradient slope depends strongly on the mass ratio of the merg- 
ing galaxies and the induced star formation during the merger. 
Typical gradients were found to be Vz = -0.22, V[p e /H] = 
-0.38 for major mergers and Vz = -0.30, V[p e /H] = —0.45 
for non-major mergers. Simulating gas-rich disk galaxy merg- 
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ers, iHopkins et al.l (j2009) also found that the amount of dissipa- 
tion, characterizing the merger-induced starbursts involved dur- 
ing the merger, is the essential factor in determining the final 
gradient slopes. Di Mat teo et al.l d2009l) investigated the shape 



of metallicity gradients resulting from dissipationless, (gas-free) 
dry mergers with the conclusion that the gradients strongly de- 
pend on the gradient shapes of the progenitor galaxies. Two 
galaxies with similar gradients will result in a flatter gradient 
in the final galaxy. However, the initial gradient can also be pre- 
served or even increased by the merger if the slope of the com- 
panion galaxy is sufficiently steep. 

The comparatively shallow metallicity gradients of our sam- 
ple (V[Fe/H] = -0.19 + 0.08) with respect to the predicted gra- 
dients from m onolithic colla pse for the most-massive galaxies 
(V z = -0.5: ICarlbergi ri984) suggest that major mergers have 
been involved in the formation of fossils, continuosly flatten- 
ing initial gradients to the present day values. Espec ially when 
compa ring our results with the simulated gradients in Kobavashi 
(2004), our gradients derived from both ULySS and lick indices 
match the distribution from major me rgers whereas no n-major 
mergers show much steeper gradients. iKobavashil (120041) argued 
that the spread in metallicity gradients and hence the formation 
of ellipticals cannot be explained by major mergers or mono- 
lithical collapse alone but both processes might be involved in 
shaping the final gradients in ellipticals. However, if fossils orig- 
inate from an at ypical, top-heavy luminosity fu nction - the failed 
group scenario (Mulchaey &Zabludoff 1999), with the magni- 
tude gap being in place a priori, i.e. the central elliptical formed 
monolithically, no L* galaxies would have been present to allow 
for subsequent major mergers after an initial monolithic collapse 
and we would expect to observe much steeper gradients. Hence 
we conclude that fossils are not failed groups but indeed formed 
through successive major mergers of L* galaxies. 
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Spola or et al. I d2009t) also found that massive galaxies with 
log cr > 2.15 show in general flatter gradients than low-mass 
galaxies, but exhibit a much larger scatter in gradient slopes. The 
flatter gradients were also attributed to galaxy mergers which 
should be the driving process in building up galaxies above this 
mass limit. The large scatter in gradient strengths was interpreted 
to be the consequence of various merging histories with differ- 
ent amounts of gas involved during the mergers. The fossil gra- 
dients in this work do not show a large spread in gradient slopes 
suggesting similar merging histories. Figures UTI and [121 show 
trends between gradient strength and central age for the Virgo 
and Fornax early-type galaxies. There also seems to be a slight 
trend between age and metallicity gradients, steeper negative 
metallicity gradients showing more positve age gradients. All 
observed fossils fall on these relations. However, because of the 
low-number statistics of our sample it is impossible to conclude 
if fossils systematicall y follow these t rends. We also compared 
our fossils to the FP of Bema rdTet al.l (l2003bl) and conclude that 
our sample follows their FP within 3<x confidence intervals, in- 
dicating that they a re dynamically relaxed sy stems, at least in 
their central parts dMendez-Abreu et al.ll2012l) . Interestingly, all 
our galaxies fall below the FP, suggesting that they are either too 
compact for given surface brightness and velocity dispersion or 
that they are too massive for a given effective radius. The latter 
seems very plausible since these galaxies are expected to consti- 
tute the mass of several U galaxies. 

5. Conclusion 

We observed spatially resolved stellar population parameters in 
a sample of fossil ellipticals. Exhibiting an average central ve- 



locity dispersion of cr = 271 + 28 km s _1 , our sample com- 
prises some of the most massive galaxies in the universe, ex- 
pected for the remnants of multiple mergers of L* galaxies. We 
find comparatively flat (V[p e /H] = -0.19 ± 0.08) metallicity gra- 
dients that suggest a hierarchical merger scenario for fossils. If 
fossils were failed groups, arising from an atypical luminosity 
function where the 2 mag gap was in place a priori, then FCGs 
would have formed monolithically, not experiencing any major 
mergers during their evolution due to the absence of L* galaxies. 
Models of monolithic collapse predict steep metallicity gradients 
(Vz = -0.5) in such a formation scenario especially for massive 
galaxies. Such steep gradients are not observed in this work. 
Hence we conclude that the observed fossils indeed formed 
through major mergers of L* galaxies. The low scatter of gra- 
dient slopes suggests a similar merging history for all galaxies 
in our sample. Our observations clearly suffer from low-number 
statistics and a much larger sample of FCGs will be needed to 
systematically quantify spatially resolved stellar population pa- 
rameters in these galaxies. Measurements of a larger sample will 
then help to probe the uniformity of the formation of these sys- 
tems. 
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